Bright cavity polariton solitons 
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The lower branch of the dispersion relation of exciton polaritons in semiconductor microcavities, 
operating in the strong-coupling regime, contains sections of both positive and negative curvature 
along one spatial direction. We show that this leads to the existence of stable one-dimensional bright 
microcavity solitons supported by the repulsive polariton nonlinearity. To achieve localization along 
the second transverse direction we propose to create a special soliton waveguide by changing the 
cavity detuning and hence the boundary of the soliton existence in such a way that the solitons are 
allowed only within the stripe of the desired width. 



The study of nonlinear and quantum effects in opti- 
cal microcavities is an area of current vigorous activity 
P S HI- I 11 particular, the strong coupling regime in 
quantum well semiconductor microcavities P, 0] leads 
to a strong and fast nonlinear response of microcav- 
ity exciton polaritons, which have been used for low 
threshold Instability [j, U IE 0] > parametric wave mixing 
@, HE, El E3, El and for the demonstration of long 
range polariton coherence supporting claims for observa- 
tion of Bose-Einstein condensation pj, G3, [ij, [IJ GJ, EH . 
Research in this area bordering condensed matter physics 
and photonics has links with many topics of current in- 
terest, which include, but not limited to, superfluidity, 
solitons, vortices, photonic crystals and quantum infor- 
mation. 

Diffraction and dispersion are important features of 
waves of any nature and the problem of their control 
stimulates significant research efforts. The recent time 
has witnessed remarkable developments into photonic 
crystals, which can be used to fully compensate diffrac- 
tion and to control its value and sign, see, e.g., [2fl [2l|. 
The ideas developed for electromagnetic waves have also 
been applied to matter waves [13], surface plasmons (23[ 
and other waves. Microcavity polaritons have been a 
part of these research efforts, with periodic potentials for 
polaritons created by the mirror patterning or by the sur- 
face acoustic waves [24], HE} . Note, that the term diffrac- 
tion is most commonly used when the beam spreading 
happens upon propagation in space. While, in the con- 
text of planar microcavities one deals with the tempo- 
ral evolution of spots of light, which is analogous to the 
spreading of quantum wave packets. This evolution is 
governed by the dispersion law linking the frequency (or 
energy) to the transverse wavenumber (or momentum). 

Diffraction and dispersion control has been compre- 
hensively discussed and used in conjunction with spatial 
and temporal soliton formation, see e.g. [26|, [U [HI . A 
particular field of these studies concerns spatial cavity 
solitons in the weak coupling regime with periodic mod- 
ulation of the detuning and/or other cavity parameters, 
see, e.g., [29j . Weak and strong coupling regime differ 
mainly in that in the latter case the dispersion curve ex- 



hibits two branches signalling mixing of light and matter 
excitations. Namely, additionally to the upper parabolic 
branch, which is largely a photonic one, the lower polari- 
ton branch appears. In our context the most relevant fea- 
ture is that this lower branch exhibits an inflection point 
where the second order dispersion changes sign [FigQJa)]. 
It is worth noting that this effect is merely evoked by the 
strong photon-exciton coupling, hence it appears even in 
a homogeneous cavity and does not require any mod- 
ulation as in the weak coupling regime. The primary 
goal of this Letter is to demonstrate how this peculiarity 
can be exploited for the formation of stable bright cav- 
ity polariton solitons (CPSs) supported by the repulsive 
exciton-exciton interaction. However, the inflection of 
the dispersion curve provides only localization in one di- 
mension and additional measures for soliton trapping in 
the other dimension have to be implemented relying on 
the soliton existence domain shift by means of the space 
dependent detuning control. 

Recently it has been suggested to use bright cavity soli- 
tons in the weak coupling regime in all optical processing 
schemes [13, HH, HH . In this context bright CPSs in the 
strong coupling regime are superior and promising can- 
didates for a practical implementation because they have 
a picosecond excitation time and require extremely low 
pump intensities of about 100 W/cm 2 , thus surpassing 
their weakly coupled counterparts by almost two orders 
of magnitude on these two crucial parameters [33]. The 
first experimental observations of dark and bright soliton- 
like structures in the strong coupling regime of semicon- 
ductor microcavities have been recently reported in [34| 
calling for more theoretical and experimental studies of 
this topic. It has been demonstrated [33] that the dark 
solitons are the only stable solutions with momenta cen- 
tered around the bottom of the lower polariton branch. 

The widely accepted dimensionless model for excitons 
strongly coupled to cavity photons is P, 0, H, HU 

d t E -i(0 2 x + d 2 y )E + [ 7c - iA - iU{y)\ E = 
= i^ + E e ik ° x , (I) 
dt^ + (70 - iA)* + i|*| 2 * = iE. 
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FIG. 1: (color online) (a) Polariton dispersion: lower (LP) 
and upper(UP) branch in the strong coupling regime, (b) 
Dispersion coefficients vs. inclination of the holding beam fco- 



Here E and \& are the averages of the photon and exci- 
ton creation or annihilation operators and polarization 
effects are disregarded. The normalization is such that 
{£lii/g)\E\ 2 and (VLr/ g)\^>\ 2 are the photon and exci- 
ton numbers per unit area. Here, fin is the Rabi fre- 
quency and g is the exciton-exciton interaction constant. 
A = (lo — w,.)/£lft describes detuning of the pump fre- 
quency u> from the identical resonance frequencies of ex- 
citons and cavity, uj r . The time t is measured in units of 
1/Q,r. 7 c and 70 are the cavity and exciton damping con- 
stants normalized to Qr. The transverse coordinates x, 
y are normalized to the value xq = c/2k z nVlR ~ 1/jto 
where c is the vacuum light velocity, n is the refractive 
index and k z = nuj/c is the wavenumber. The normalized 
amplitude of the external pump Eq is related to the phys- 
ical incident intensity Ii nc as |-Bo| 2 = gicIinc/h-^^R 2 
(l3| . fco is the transverse wavenumber of the pump, with 
fco = corresponding to normal incidence. As a guide- 
line for realistic estimates one can use parameters of a 
microcavity with a single InGaAs/GaAs quantum well: 
Kl R ~ 2.5 meV, hg ~ 10" 4 eV/^m 2 , see 0, @|. As- 
suming the relaxation times of the photonic and exci- 
tonic fields to be 2.5 ps gives 7 C) o — 0.1. In accordance 
with this set of parameters the normalized driving ampli- 
tude |-Eo| 2 = 1 corresponds to an external pump inten- 
sity of about 10 kW/cm 2 . Optical bistability appears for 
\E P \ ~ 0.1 or about 10 2 -10 3 W/cm 2 . U(y) is the spatially 
varying shift of the cavity resonance which represents a 
potential for photons. The spatial inhomogeneities of the 
cavity and quantum wells are assumed negligible, so that 
long range spatial effects can be observed. As, for exam- 
ple, in the recent demonstration of the 50/im in diameter 
anti-phase polariton 7r-state in a microcavity with peri- 
odically patterned mirrors (24|. 

To understand the physics of CPS formation and to 
identify domains of their existence we first expand the 
linear polariton dispersion around the pump wavevector, 
and arrive after an inverse Fourier transformation at a 
single equation with an effective nonlinearity and multi- 
ple dispersion terms. We assume U — and look for the 
linear polariton eigenstates, i.e., the photons trapped in 
the pump free cavity (Eq = 0, A — > — u) r /Q,n) coupled 
to noninteracting excitons. Seeking solutions of Eqs. (1) 
in the form E,^ = ek,ipke~ tnt+tkx gives the eigenvalue 



problem: 

»n ± (fc)A = ( lp + 7c + - l la ) p k , 

(2) 

where pk — {ek,ipk} is the polariton basis vector and 
n±(k) are the eigenfrequencies, corresponding to the 
upper- (+) and lower- (-) branch of the polariton disper- 
sion curve, respectively, see Fig. [TJ The inflection point 
<9 2 fi_ = is located at k — kd, where second order dis- 
persion disappears. 

The standard next step here would be to disregard the 
upper branch and to consider the fco values close to zero, 
see, e.g. [H, H, [E HI- I n order to embrace the inflec- 
tion point and to account for the varying dispersion we 
modify this technique and include higher-order dispersion 
terms. To achieve this we perform the Fourier transform 
{E, ~ J a(t, k)pke tkx dk, where a(t, k) is the Fourier 
amplitude of the fc-th component and pk corresponds to 
the lower polariton branch. We assume that the spec- 
trum of the polariton wavepacket is centered around fco 
and expand f2_(fc) up to the fourth order in \k — ko\. The 
slowly varying amplitude of the polariton wavepacket is 
then A(t,x) = a(t, k)e l(k - k °^ x d(k - fc ) and obeys: 

id t A + iD x d x A + D 2 d 2 x A - iD a d^A - D^A- 
-6A-Z\A\ 2 A = ir,E , (3) 

where 6 = f2_(fco) — uj/Qr, Red is the polariton 
detuning from the pump frequency and Im S is the 
loss. D\ = dk$l-\k is the transverse group velocity 
and D 2 = (l/2)9 2 r!_U , D 3 = (1/6)S|0_ | feo , D 4 = 
(l/24)(9 4 f2_|fc are the dispersion coefficients [Fig. [TJb)]. 
£ = Hkol 4 / (\ek„\ 2 + \ipk Q \ 2 ) is the effective nonlinearity 
and V = e k /(\ek \ 2 + \^k )\ 2 )- For fc varying from zero 
and towards, (but not too close to kd) D% > dominates 
over D m> 2 [Fig. [ljb)]. In this case Eq. ((3]) coincides 
with the well-known weak-coupling Lugiato-Lefever (LL) 
model |35j |. 

In the LL models applied for the lower polaritons the 
dispersion is always normal, D 2 > 0, and the nonlinearity 
is repulsive, therefore the only stable CPSs are the dark 
ones [331 ] . while bright solitons exist, but are unstable 
[Fig. [2{a,b)]. However, D 2 (k) crosses zero at fco = kd and 
the Z?3-term becomes the leading one around this point. 
Then the I?3-term drops and the polariton dispersion is 
determined by the competing D 2 and D4 terms. Both 
D 2 < and D4 > favor the existence of bright solitons 
for the repulsive nonlinearity and therefore we expect 
to find them for fco > kd- Moreover, it is evident that 
because of D\ ^ these solitons are forced to move. 

While the above considerations are useful for qualita- 
tive understanding, we return to the original model (1) 
with d y — in order to compute ID bright CPSs. In 
seeking for soliton solutions we have to take their motion 
into account: E(t,x) = E(£)e ik ° x , *(*,») = S&(£)e ikoX , 
£ = x — vt, where v is the velocity (yet to be determined) 
close, but not equal, to D\, see Fig. 2(a). E and ^> 
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FIG. 2: (color online) Maxima of the excitonic components of 
the solitons vs + (left column) and the soliton profiles (right 
columns) for different fco- Black lines in the left column mark 
the peak amplitudes of the spatially homogeneous solutions 
(HS) and the red ones of the solitons. The full/dashed lines 
in the left column correspond to stable/unstable solutions re- 
spectively. Eo and A values are: (a,b) 0.175 and —0.7; (c,d) 
0.149 and -0.39; (e,f) 0.19 and -0.05 . 

obey (2fc - v)d 6 E - id\E + ( 7c - iA + ik$)E = W + E 

and -vd^ + ( 7o - iA)* + i|*| 2 * = iE. Bright CPS 
solutions of the above model are found by a modifica- 
tion of the Newton method allowing to treat v as an 
unknown variable. The maximum of v occurs at k = kd, 
then the soliton velocity decreases with fco increasing 
and approaches zero, see Fig. 3(a). This is in remark- 
able contrast to what happens if dispersion is parabolic, 
where the soliton velocity continuously increases with fco. 
For our choice of parameters v must be multiplied by 
4 x 10 6 m/s to give the physical velocity. It implies that 
the CPS with v = 0.25 will traverse across the typical 
distance of 100/im in lOOps. This is 40 times larger than 
the polariton life time. Hence, once excited, the solitons 
have enough time to shape and be observed during their 
motion. 

Plots illustrating the dependence of the maximal soli- 
ton amplitude from the pump field and the associated 
bistability curves for the homogeneous (d XjV = 0) solu- 
tions are shown in the left column of Fig. 2. The right 
column shows the corresponding transverse soliton pro- 
files. For fc = we are reproducing the results of [33| 
with bright CPSs being unstable. For fco just below and 
above kd the branch of bright CPSs starts bending and 
islands of their stability start to emerge and expand with 
fco increasing. 

To reveal the influence of different dispersion orders on 
the CPSs we plotted the dependence of the CPS width 
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FIG. 3: (color online) (a) Diamonds show numerically com- 
puted soliton velocity vs. fco- Dashed line shows Ui(fco). (b) 
Diamonds show the soliton width vs fco- Dashed lines show 
the appropriately scaled 7 \/D m , where m = 2, 3, 4. 




FIG. 4: (color online) (a) Soliton moving along the potential 
attractive for linear polaritons: U = 6exp[— (j//2) 8 ], b = 2, 
Eo — 0.18. Other parameters as in Fig. 2(e,f). (b) Soliton 
moving along the repulsive potential b — —1, Eo = 0.213. 



on the inclination parameter fco in Fig. [3jb) . Since not 
only the dispersion coefficients vary with fco, but also 
Re 6 (which strongly influences the CPS width), we have 
kept the latter fixed at Re 5 = —0.3, by adjusting the 
detuning A respectively. For large fco the excitonic part 
of CPSs dominates over the photonic component, see Fig. 
2. That is why the CPSs can be much narrower than the 
pure photonic cavity solitons in the weak coupling regime 
and may attain widths well below the ones allowed by the 
photonic dispersion, although this happens at the price of 
oscillatory tails. The CPS width is expected to scale with 
the dominant dispersion coefficients as ~ \/D m (m = 
2, 3, 4). Fig. [3^b) compares the numerically found soliton 
width (diamonds) with the scaling given by the different 
dispersion orders (dashed lines). Third-order dispersion 
describes well the soliton width for relatively small fco 
where D2 — 0. Further inclination brings D2 on the 
top, while the third-order dispersion vanishes. A further 
increase of fco brings fourth-order dispersion into play, 
which starts to compete with £>2- 

In order to use CPSs as information bits, it is desir- 
able to have them robustly localized in both transverse 
dimensions. However in the homogeneous cavity this is 
impossible to achieve, since the signs of polariton disper- 
sion are opposite in orthogonal directions. To solve this 
problem we propose to shift the cavity detuning in such 
a way that outside the waveguide area the soliton exis- 
tence conditions are not satisfied. This is particularly 
easy to achieve if the operating frequency is not far from 
the boundary of the existence domain. With this ap- 
proach, it does not matter if one uses U(y), which traps 
linear waves, or U(y), which repels them. Thus using 
of strongly nonlinear polariton waves extends the list of 
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FIG. 5: (color online) (a) Excitonic part of the soliton spec- 
trum (IV'fcl) of the the soliton shown in Fig. 2(d). (b) The 
density plot is the 2D generalization of \"4>k\ from (a). The 
slope of this plot shows dependence of the soliton frequency 
on k. The full line is LP dispersion. 

the possible trapping techni que s already explored for lin- 
ear cavity polariton [l?], HH ■ Fig. QJa) demonstrates 
soliton guiding in the trapping potential and Fig. |3)Jb) 
in the repelling potential. The trapping potential looks 
preferential for the soliton guidance, because it reduces 
soliton interaction with the quasi- linear waves. The lat- 
ter are far more pronounced in the case of the repelling 
potential, see Fig. [3Jb). Practical excitation of CPSs is 
no different from their weakly coupled counterparts and 
requires short term application of the addressing beam 
in addition to the pump Eq [U . 



In conclusion, we have shown that stable moving ID 
cavity polariton solitons exist providing the lower polari- 
tons are used near or beyond the inflection point of their 
dispersion. We have also demonstrated the continuous 
transformation between cavity soliton polaritons shaped 
by the different dispersion orders, including the regimes 
where the influence of the usual second order dispersion 
is negligible. We have proposed a mechanism for their 
guiding and localization in the second transverse dimen- 
sion. Our results form the basis for future studies of 
soliton-polariton logic and processing schemes. 

Note. During the revision of this manuscript the re- 
lated experimental observations of moving quasi-localised 
polaritons have been published [37j . Brief preliminary 
comparison of dispersion plots from (37j with our re- 
sults is very encouraging. To reveal the soliton fre- 
quency content we take the Fourier transform of its ex- 
citonic component: tp^ = J x ty(x — vt)e lk ° x+lkx dx = 

f(k - k )e M ^- k \ where f(k) = J^V^e^d*. \^ k \ 
is plotted in Fig. 5(a). Apart from the infinite peak 
at the at fco corresponding to the background, the soli- 
ton spectrum proper shows 3-peak structure, such that 
2k p = ki + k s , where k p « ko k ( i- The idler field ki 
was externally seeded and the signal k s was generated in 
Ref. [I?]], while both of them are the integral parts of 
our soliton solution. Soliton frequency calculated from 
ijjk is a straight line shifted upwards from the dispersion 
of linear polaritons (Fig. 5(b)). Figs. 5(a,b) essentially 
match the corresponding data from [37| . 
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